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Abstract 

The propagation of small-amplitude inhomogeneous plane waves in 
an isotropic homogeneous incompressible Mooney-Rivlin material is 
considered when the material is maintained in a state of finite static 
homogeneous deformation. Disturbances of complex exponential type 
are sought and all propagating inhomogeneous solutions to the equa- 
tions of motion are given, as well as the conditions for linear, ellip- 
tical, or circular polarization. It is seen that a great variety of solu- 
tions arises. These include some original solutions, such as circularly- 
polarized plane waves which propagate with an arbitrary complex 
scalar slowness, or linearly-polarized waves for which the direction of 
propagation is not necessarily orthogonal to the direction of attenua- 
tion. Throughout the paper, geometrical interpretations and explicit 
examples are presented. 
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1 Introduction 



The theoretical study of elastic plane waves has generated a large litera- 
ture in the area of finite elasticity. One type of plane waves which is of 
practical interest is that of motions propagating in a finitely and homoge- 
neously deformed elastic material [T], because the predeformation can be 
used to model, at least locally, the anisotropy of many mechanical and ge- 
ological structures. A review of these topics can be found in a textbook 
by Ie§an [2]. Once it is homogeneously deformed, an elastic body presents 
three privileged orthogonal directions, namely those of the principal axes of 
the static deformation. Waves propagating in one of these directions are 
called 'principal waves' |3] and their properties are well known and relatively 
easy to establish, because the principal directions offer a natural rectangu- 
lar Cartesian coordinate system. Homogeneous plane waves propagating in 
an incompressible elastic material are necessarily transverse (the direction of 
propagation is orthogonal to the plane of polarization). Therefore, for such 
waves propagating in a nonprincipal direction, another rectangular Cartesian 
coordinate system proves to be useful, that formed by the direction of prop- 
agation and by two orthogonal directions in the plane of polarization. Currie 
and Hayes |1] proved that the Mooney-Rivlin form for the strain energy den- 
sity (which is used to model the mechanical behaviour of rubber [HE]), is 
the most general one for which homogeneous plane waves, be they of finite 
or of small amplitude, may propagate in any direction for an arbitrary finite 
static homogeneous predeformation. Later, Boulanger and Hayes [3 [8] stud- 
ied in great detail finite-amplitude homogeneous plane waves in a deformed 
Mooney-Rivlin material. 

Homogeneous plane waves are such that the planes of constant phase 
(orthogonal to the direction of propagation) are parallel to the planes of 
constant amplitude (orthogonal to the direction of eventual attenuation). 
However, for a variety of physical problems, an attenuation of the amplitude 
occurs in a direction distinct from the direction of propagation. In those 
cases, a combination of 'inhomogeneous' plane waves is introduced, usually 
in the form e""^^ ■''cosu;(S' • x — t)a, where u is the frequency, and S', S", 
a are vectors. This procedure was successfully applied to various interfacial 
problems, such as reflected and refracted waves, Rayleigh waves. Love waves, 
Stoneley waves, Scholte waves, etc. When the directions of propagation 
(that of the vector S'), of exponential attenuation (that of the vector S"), 
and of polarization (that of the vector a) are orthogonal with respect to one 
another, they form the basis for a rectangular Cartesian coordinate system 
in which the incremental equations of motion can be written and solved for 
small-amplitude [9l [10] as well as for finite-amplitude [TT] inhomogeneous 



2 



waves in a deformed Mooney-Rivlin material. When these directions are not 
orthogonal, the equations of motion become much harder to solve. This is 
where the algebra of complex vectors, also known as 'bivectors', can play 
a very useful role. The use of bivectors in order to describe inhomogeneous 
plane waves is made clear in a textbook by Boulanger and Hayes [12] . Indeed 
bivectors, that is vectors with a real part and an imaginary part, can be used 
to describe the polarization of a wave through the 'polarization bivector' A = 
A' + iA", as well as its propagation and attenuation through the 'slowness 
bivector' S = S' + iS". Hence an inhomogeneous plane wave is modelled as 
being proportional to the real part of the expression Ae*'^*^^ ''"*-'. 

Previously, the propagation of small-amplitude inhomogeneous plane waves 
in a Mooney-Rivlin material subjected to a finite static biaxial homogeneous 
deformation has been considered first by Belward [13], and later by Boulanger 
and Hayes [H] , using bivectors. In this paper, we consider the propagation of 
infinitesimal inhomogeneous plane waves in a Mooney-Rivlin material which 
is maintained in a state of arbitrary triaxial finite static homogeneous de- 
formation. Most results obtained here are a generalization and an extension 
to the case of inhomogeneous waves of results established by Boulanger and 
Hayes [7] for the propagation of finite-amplitude homogeneous plane waves in 
a homogeneously deformed Mooney-Rivlin material. However, this extension 
is only possible when the amplitude of the wave is considered small enough 
to allow linearization. This restriction is due to the fact that finite-amplitude 
inhomogeneous plane waves can propagate in an incompressible elastic ma- 
terial only when they are linearly-polarized [15]. For infinitesimal inhomo- 
geneous plane waves, no such restriction applies and elliptical polarization 
is possible. The purpose of this paper is to find all inhomogeneous small- 
amplitude plane waves of complex exponential type travelling in a deformed 
Mooney-Rivlin material, and to establish the conditions for linear, elliptical, 
and circular polarization. Within this context, a much greater number of 
solutions are found for inhomogeneous waves than for homogeneous waves. 
For instance, elliptical polarization is possible for homogeneous waves only in 
two special directions (the 'acoustic axes' [7]), whereas it is possible in other 
directions for inhomogeneous waves. Also, for certain inhomogeneous waves, 
the 'complex scalar slowness', which is the counterpart of the inverse of the 
speed for homogeneous waves, may be arbitrarily prescribed. In general, for 
a given orientation of the plane containing the directions of propagation and 
of attenuation, there is an infinity of inhomogeneous wave solutions. 

The paper is organized as follows. First (section we recall the basic 
equations describing the Mooney-Rivlin material. Then (section 12]) we write 
the equations governing the motion of a small-amplitude disturbance in a ho- 
mogeneously deformed Mooney-Rivlin material. The incremental equations 
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of motion and the incompressibility constraint are given. In order to solve 
these equations, we seek solutions of complex exponential type, which are 
presented in section HI The slowness bivector S is introduced; when the real 
and the imaginary parts of S are not parallel, the wave is inhomogeneous, be- 
cause the directions of propagation and of attenuation do not coincide. The 
amplitude bivector A is also introduced; depending on whether or not the real 
and the imaginary parts of A are parallel, the wave is linearly or elliptically 
polarized, respectively. The special case of circular polarization corresponds 
to the 'isotropy' of A, that is A ■ A = [I6]. Using these quantities, we then 
derive the incremental equations of motion for inhomogeneous plane waves 
of complex exponential type propagating in the deformed Mooney-Rivlin 
material. 

Next (section^ and E]) the propagation of such disturbances is investi- 
gated. Different sub-cases arise, depending on whether or not the slowness 
complex vector S is isotropic, and on whether the waves are circularly, ellip- 
tically, or linearly polarized. For each case, the general solution is provided 
(amplitude and slowness bivectors, incremental pressure, wave speed) as well 
as various explicit solutions to the incremental equations of motion. 

2 Basic equations 

The Mooney-Rivlin material is a homogeneous isotropic hyperelastic solid, 
for which the strain energy E per unit volume is given by ^ 



where C, D are material constants, and /, // are the first and second in- 
variants of the left Cauchy-Green strain tensor B. In order to satisfy the 
strong ellipticity condition, the constants C, D are such that [T71 18] either 
C > 0, D > or C > 0, D > 0. When D = 0, the material is said to be 
'neo-Hookean', but this possibility is not considered in this paper. 

The strain tensor B is related to the deformation gradient F through 
B = FF""", and its first two invariants / and // are defined by / = trB and 
211 = (trB)2 -tr(B2). 

Because the Mooney-Rivlin material is incompressible, we have at all 
times. 



Finally, the Mooney-Rivlin constitutive equation is derived from (12. ip 



2S = C(/-3) + D(//-3), 





(2.2) 



and (O as ^ 



cr = — 



{p - DII)1 + CM - DM 



1 



(2.3) 
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where cr is the Cauchy stress and —pi a hydrostatic pressure, to be deter- 
mined from the equations of motion, and initial and boundary conditions. 

Throughout the paper, we assume that no body forces are apphed, so 
that the equations of motion reduce to 

diver = (2-4) 

Here, x is the position at time t of a point of the deformed body which was 
at X in the undeformed state, and p is the constant mass density. 

3 Small motions superposed on a large ho- 
mogeneous deformation 

In this Section we give the incremental equations of motion corresponding 
to the superposition of a small-amplitude motion upon a large homogeneous 
strain. First, we assume that the Mooney-Rivlin material is subjected to a 
finite pure homogeneous deformation, for which a point initially at X in the 
rectangular Cartesian coordinate system (O, i, j, k) moves to a point at x in 
the same system, with extension ratios Aq (a = 1, 2, 3) and principal axes of 
deformation along the axes of the coordinate system, so that 

Xa = ^aXa (a = 1, 2, 3; no sum). (3.1) 

The extension ratios are assumed to be distinct, constant, and without 
loss of generality, to be ordered as Ai > A2 > A3. 

The deformation gradient F and strain tensor B corresponding to this 
finite static deformation are constant and given by F = diag (Ai, A2, A3) and 
B = diag (A?, A^, A^). Then the tensor T, defined by P [7] 

T„/3 = 0, a 7^ /3, 

Taa = -V + D{\f + A^2 + \f) + C\-^ - DX-\ (a = 1, 2, 3; no sum) 

(3.2) 

where p is constant, is the constant Cauchy stress required to sustain the 
static deformation f l3.ip . 

Boulanger and Hayes [7] have shown that two axes, namely the 'acoustic' 
axes of the B~^-ellipsoid, exhibit remarkable properties with respect to the 
propagation of finite-amplitude homogeneous plane waves in such a deformed 
Mooney-Rivlin material. The acoustic axes are the only directions in which 
circularly-polarized finite-amplitude homogeneous plane waves may propa- 
gate. Their directions n='= are defined in terms of the stretch ratios Ai, A2, A3, 
and of the unit vectors i and k, which are in the directions of the principal 
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axes of B (or B ^) corresponding to Ai and A3, respectively. These directions 
are given by 



= cos 0i ± sin 0k, where cos = \ — -39 , sin — '^^ 




Ag^ - A^^ y Ag^ - A^^ 

(3.3) 

and are independent of the material parameters C and D. They lie along 
the normals to the planes of central circular sections of the B~^-ellipsoid fl2\ 
§5.7]. Their directions also appear in the Hamilton cyclic decomposition 
of the B~^ tensor as, 

B"^ = A^^l - ^{Xf - Ar^)(n+ ® n- + n- ® n+). (3.4) 

Now we consider a further deformation, possibly time-dependent, in which 
the particle at x moves to x such that 

x = x + eu(x, t). (3.5) 

Here the displacement eu is a vector depending on x and on the time t, and 
e is a small parameter. Throughout the paper, we neglect terms of second 
and higher order in e. 

Now if the deformation gradient F, the strain tensor B, the invariants J, 
//, the pressure p, and the stress tensor T corresponding to the deformation 
(13.51) are expanded in powers of e around their their value in the static state 
of homogeneous deformation, we get 

F = F + eF* + . . . , 1 = B + eB* + . . . , 

/ = / + £/* + ..., 1! = 11 + eir + (3.6) 

p = p + ep* + T = T + eT* + . . . 



Then the equations of motion (12. 4p written for T and x take the form [18] 

divT = p{d^u/dt^), (3.7) 

when we retain terms up to order e. Explicitly, for the Mooney-Rivlin case, 
they are found to be (see also Belward jj^, Hayes and Horgan [19] for the 
biaxial case): 

-p*, + CXjui^u + {CXI + DK')U1,22 + (CXI + DA^^)mi,33 

-DA^%,12 - -DA3"^U3,i3 = pill, 

-p% + {CXj + DX2^)U2,U + CXIU2,22 + {CXj + DX^^)U2,33 

-DX^^Ui^U - DX^^U3^23 = PU2, 

-p% + {CXI + DXf)u3,n + {CXI + DXf)u3,22 + CXlu3,33 

-DX^'^Ui^l3 - -DA^^M2,23 = PU3, 

6 



(3.8) 



where commas and dots denote differentiation with respect to position x and 
time t, respectively. Thus, for example, Ui^2 = dui/dy. 
Finally, the incompressibility constraint (12. 2 p yields [18] 

+ M2,2 + M3,3 = 0. (3.9) 

4 Vibrations of complex exponential type 

In order to solve the equations governing the incremental motion, we assume 
that the displacement eu and the incremental pressure ep* are of the form 

u = |{Ae*"(^-^-*) + c.c.}, p* = |{zwQe*^(^''~*) + c.c.}, (4.1) 

where A = A' + iA" is a complex vector, called the 'amplitude bivector' |12] . 
uj is the real frequency, Q is a scalar, S = S' + iS" is the 'slowness bivector', 
and 'c.c' stands for 'complex conjugate'. 

An ellipse may be associated with a bivector, as the ellipse for which the 
vector corresponding to the real part and the vector corresponding to the 
imaginary part are conjugate vectors. That is, to a bivector D = D' + iD" 
say, we associate the ellipse described by the point M such that OM = 
D' cos e + D" sin ^, < ^ < 2tt. 

When the real and imaginary parts of A are parallel, that is, when the 
ellipse of A degenerates into a segment, the wave is linearly-polarized along 
their common direction. Otherwise, the wave is elliptically-polarized, and 
the ellipse of polarization is the ellipse of A. The special case of circular 
polarization corresponds to the 'isotropy' of A [Tn| 

A-A = 0. (4.2) 

The planes defined by S' ■ x = constant are called planes of constant 
phase, and the planes defined by S" • x = constant are the planes of constant 
amplitude. When the real part S' and the imaginary part S" of S are not 
parallel, the wave is said to be inhomogeneous. In that case, we introduce 
the 'directional ellipse' of the slowness bivector S, defined as follows. Let 
m and n be unit vectors along the respective major and minor semi-axes 
of the ellipse of S, and m be the aspect ratio of the ellipse of S. Then the 
directional ellipse of the slowness bivector S is the ellipse of the 'propagation 
bivector' C, defined by 

C = mm + in, with m-m = n-n = l, m-n = 0. (4.3) 

Then S may be written as [16] S = A^C = N{mm + in), where = A^' + 
iN" is a complex number, called the 'complex scalar slowness'. Then the 
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directions of propagation and attenuation are the directions of S' and S', 
respectively, given by 



S' = mN'm - N"n, S" = mN"m + N'n. (4.4) 

Note that 

S' . S" = {m^ - 1)N'N", (4.5) 

and therefore, the planes of constant phase are orthogonal to the planes of 
constant amplitude either when m = 1 or when is purely real or purely 
imaginary. 

So, we may write solutions of the form ( 14.11) in terms of C as 

u = i{Ae^'^(^^-''-*) + c.c.}, p* = |{iwArPe^'^(^*^''-*) + c.c.}, (4.6) 

where P = N~^Q. Now we substitute these expressions into (13. 8p and ( 13. 9p . 
First, equation (13. 9p imposes the condition 

A • C = 0. (4.7) 

A geometrical interpretation [121 §2.4] of this equation is that the orthogonal 
projection of the ellipse of the bivector C onto the plane of the bivector A is 
similar and similarly situated to the ellipse of A, rotated through a quadrant. 
The amplitude and propagation bivectors are said to be 'orthogonal'. For 
homogeneous plane waves, the equation reduces to A ■ n = (where n is 
a real vector in the direction of propagation), which simply means that the 
polarization of the wave is transverse. 

Next, the equations of motion ( 13. 8 p yield 

- PC + C(C-BC)A + D[{C ■ C)M-^A - (A-B~^C)C] = piV'^A. (4.8) 

Taking the dot product of this last equation with C, and using equation 
(grp . yields 

-P(C-C) = 0. (4.9) 

In conclusion, the propagation of small-amplitude waves of complex ex- 
ponential type in a homogeneously deformed Mooney-Rivlin material is gov- 
erned by the following equations, 

-PC + C(C-BC) A + D[(C ■ C)l - C ® C]B-i A = pN'^A, 
A-C = 0, P(C-C) = 0. ^ ^ 

Now, we treat in turn the case where C is isotropic (C ■ C = 0) and the 
case where C is not isotropic (C ■ C 7^ 0). 
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5 Propagating evanescent solutions, C • C = 



Here it is seen that corresponding to any isotropic bivector C, there exists an 
infinity of hnearly, circularly, and elliptically polarized inhomogeneous plane 
wave solutions. 

First, we note that when C • C = 0, we have m = 1 in (14. 3p and we 
deduce from fl4.5l) that the planes of constant phase are orthogonal to the 
planes of constant amplitude. 

Next, the equations (14. 8 p and (14. 9 p reduce to 

-PC + C(C-BC)A - D{A-M'^C)C = pN~^A, 

A -0 = 0, C-C = 0. ^ ' 

Equations (15. IL n allow us to decompose A as [121 §2.9] 

A = aiC + a2CAC, (5.2) 

where ai, a2 are real constants and C is the complex conjugate of C. Note 
that because C is isotropic, it is written as C = m + in ( (I4.3p with m = 1) 
and so, C A C = 2m A rii is parallel to a real vector. Substituting (15. 2p in 
(15.11) ^ leads to two different types of solutions. 



(i) First type of solution: linearly and elliptically polarized waves 

Here the amphtude bivector A is given by A = aiC + a;2C A C, 02 7^ 0, 
where ai, are real constants. The corresponding complex scalar slowness 
is given by 

= C(C-BC), (5.3) 
and the incremental pressure is given by (14. 6p . where 

P = -aiD(C-B~^C) - a2D{C A C) ■ B-^C. (5.4) 

Note that in this case, A ■ A = — q;|(C • C)^ 7^ 0, and thus the wave is 
not circularly-polarized. When ai = 0, the wave is linearly-polarized in the 
direction of CAC. When ai 7^ 0, the wave is elliptically-polarized. Also, note 
that (15.30 can be written in terms of the slowness bivector Sasp = CS-BS. 
The imaginary part of this equation yields S' ■ BS" = 0, which means that 
the direction of the normal to the planes of equal phase and the direction 
of the normal to the planes of equal amplitude are conjugate directions with 
respect to the B-ellipsoid. This condition has been established previously for 
linearly-polarized finite-amplitude inhomogeneous plane waves propagating 
in a deformed Mooney-Rivlin material with an isotropic slowness bivector 
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(ii) Second type of solution: circularly-polarized waves 

Here the amplitude bivector A is given by A = aiC, ai 7^ 0, where ai is a 
real constant, and the incremental pressure is given by f l4.6p . where 

P = ai[C(C-BC) - D(C-B-^C) - pN-% (5.5) 

For this wave, the complex scalar slowness N is arbitrary (a similar situation 
was encountered by Boulanger and Hayes [H]); thus the displacement, which 
is the real part of e{aiC exp iuj{NC ■ x — t)} is independent of the material 
constants C and D. Also, A • A = a^(C ■ C) = and the wave is circularly- 
polarized. 

Because there is an infinity of choices for an isotropic bivector C, there 
is a triple infinity of propagation, attenuation, and polarization directions 
for linearly (type (i), ai = 0), elliptically (type (i), ai 7^ 0), and circularly 
(type (ii)) polarized inhomogeneous plane waves of complex exponential type, 
provided the slowness bivector is isotropic. This is in sharp contrast to 
elliptically and circularly polarized homogeneous plane waves, which can only 
propagate along an acoustic axis [7]. 

Now we present explicit examples of the two types of solution correspond- 
ing to an isotropic propagation bivector C. 

Example 1: Waves with an isotropic slowness bivector 

We choose an isotropic bivector C and write the corresponding two types 
of solutions. One set of solutions consists of elliptically-polarized waves, the 
other set consists of circularly-polarized waves travelling with an arbitrary 
speed. Let C = (3i + 5zj + 4k)/5. 

For waves of type (i), we have, using f l5.2p . A = (3a;i + 4^0:2)1 + 5iaij + 
(4ai - 3ia2)k, pN-^ = C(9A? - 25\l + 16A2)/25 and P = -D[{9\j - 25A^ + 
16Ai)ai + 12i{X^^ - Xf)a2]/5. 

If the primary finite deformation is such that 9A^ — 25A2 + 16A| > 0, 
then N is real and we can write the following solution to the incremental 
equations of motion in a deformed Mooney-Rivlin material f l3.8p : 

Ml = e-'^l^l^{3ai cosu;[|A^|(3x + Az)/5 - t] - 4^2 sinw[|A^|(3x + 4z)/5 - t]}, 
U2 = -5aie-'^l^l^sinw[|A^|(3a; + 4z)/5 -t], 

M3 = e-'^l^l^{4ai cosa;[|A^|(3a; + 4z)/5 - t] + 3^2 smuj[\N\{3x + 4:z)/5 - t]}, 
p* = lDuj\N\e-'^\^\y{ai{9Xf-25X2^ + IQXf) smuj[\N\{3x + Az)/5 - t] 

+12a2(A^2 - Ag^) cosa;[|A^|(3x + 4z)/5 - t]}, 

(5.6) 
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where cj, ai, are arbitrary (0^2 7^ 0), and 

\N\ = 5y/p/[C(9A2-25Ai + 16Ai)]. (5.7) 

This wave travels in the direction of 3i + 4k and is attenuated in the y- 
direction. If 9A^ — 25\l + I6A3 > 0, then a solution is given by 

ui = e-'^l^l(3^+^^)/5{3aicosu;[-|A^|y-t] - 4^2 sinu;[-|A^|?/ - t]}, 
U2 = -5aie-^l^l(3^'+^")/5 smu[~\N\y - t], 

Us = e-'^l^l(3^+^^)/5{4aicosa;[-|A^|?/-t] +3a2 smuj[-\N\y - t]}, 
p* = |Du;|Ar|e-"l^l(=^"+^")/^{ai(9Ar' - 25Xf + IQXf) smuj[-\N\y - t] 

+12a2(A^2 _ ^-2) cosuj[-\N\y - t]}, 

(5.8) 

where u, ai, 0^2 are arbitrary (0:2 7^ 0), and |A^| is again given by (15. 7p . This 
wave travels in the direction of — j and is attenuated in the direction of 3i+4k. 
In both cases (15. 6 p and (15. 8p . the wave propagates with speed |A^|~^ where 
|A^| is given by (15. 7p . and is elliptically-polarized (two conjugate radii of the 
ellipse are ai(3i + 4k) and 5Q;ij + 0:2 (4i — 3k)). 

For waves of type (ii), we have: A = ai(3i + 5ij + 4k), pN~'^ is arbitrary 
and P = ai[C{9Xl - 25X1 + 16Ai) - D{9Xf - 25X2^ + IGA^^) - pN-^]/5. 
For real, the corresponding solution is given by 

cosaj[A^(3x + 4z)/5 — t], 
U2 = -5aie-'^^^sinu;[A^(3a; + 4^)/5 -t], 

cosa;[iV(3x + 4z)/5 — t], 
p* = laiuNe"'^y{pN-'' - C{9Xl - 25X1 + 16A^) 

+D{9X^^ - 25X2^ + l6Xf)} sina;[A^(3a; + 4z)/5 - t], 

(5.9) 

where 0:1,0; and are arbitrary. This wave travels in the direction of 3i + 
4k with an arbitrary speed A^~^, is attenuated in the ^/-direction, and is 
circularly-polarized (two orthogonal radii of the circle are a;i(3i + 4k) and 
5aij). The displacement field eu does not depend on the constants C and 
D. It is easily checked that the solutions (15. 6p . (15. Sp . and (15. 9p satisfy the 
general equations of motion (13.80 . 

6 Propagating evanescent solutions, C • C 7^ 

In this Section, we consider the case where the bivector C is not isotropic. 
A great variety of solutions is uncovered, and a systematic method of con- 
struction and classification for linearly, circularly, and elliptically polarized 
waves is presented. 
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When C ■ C 7^ 0, dHODs implies that P = and therefore p* = 0. 
Equation (14. 101) 1 then reduces to 



{C(C-1C)1 + D[{C ■ C)l - C ® C]B-i}A = pN-'A. (6.1) 

Following [7], we symmetrize this equation, using (14 .Tp . to give the equivalent 
form 

n[C(C-BC)l + D(C-C)B-^]nA = pAr-2A and A ■ C = 0, (6.2) 
where we have introduced the 'complex' projection operator 11, defined by 

This operator generalizes the 'real' projection 1 — n ® n upon the plane 
n-x = [7J, and has the following properties: 11^ = 11, IIC = and 

nA = A. 

By inspection of (16. 2p . we see that solving the equations of motion, once 
C is prescribed, is equivalent to finding the eigenbivectors A of the tensor 
n[C(C-BC)l + D{C ■ C)B-i]n such that A ■ C = 0, and their correspond- 
ing eigenvalues pN'"^ . This procedure is analogous to that used in [Z] for 
finite-amplitude homogeneous plane waves, with the replacement of their 
real vectors of propagation n and polarization a, and their speed c with the 
bivectors C and A, and the quantity A^~^, respectively. 

First of all, we compute the possible eigenvalues pN~'^. 



6.1 Secular equation 

The equations (16. 2p admit solutions, provided that 

det(n[C(C-BC)l + D{C ■ C)B-^]n - pA^^^l) = 0. (6.4) 

Equation (16. 4p is the classical secular equation for inhomogeneous plane 
waves. Because det 11 = 0, then pN~'^ = is one root of (16. 4p . The two 
other eigenvalues pN^^ (say) are the roots of the quadratic 

[pN-^ - C(C-BC)]2 - D[{C ■ C)trB-i - (C-B-^C)] [piV-^ _ C(C-BC)] 

+ D2(c.C)(C-BC) =(6.S) 

Explicitly, pN^"^ are given in terms of C alone as 

2piV±2 = 2C(C-BC) + D[{C ■ C)trB-i - (C-B-^C)] 

± D^[{C ■ C)trB-i - (C-B-iC)]2 - 4(C ■ C)(C-1C). (6.6) 
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These last quantities are the two eigenvalues corresponding to the amplitude 
bivectors which are orthogonal to C. Equations fl6.6p generalize the corre- 
sponding equations written for homogeneous waves by Boulanger and Hayes 
[8], by replacing their vector of propagation n and speed c by C and A^~^, 
respectively. Now we establish in turn the conditions for linear, circular, and 
elliptical polarization. 

6.2 Linearly-polarized waves, C • C 7^ 

Here we prove that linearly-polarized inhomogeneous waves with a non- 
isotropic slowness bivector can propagate in a deformed Mooney-Rivlin ma- 
terial only when they are polarized in a principal direction. 

When the amplitude bivector A is parallel to a real unit vector a (linear 
polarization) , the incompressibility constraint (14. 7p written as a ■ S = im- 
plies that the directions of propagation (S') and attenuation (S") are both 
orthogonal to a. Hence, we write a = aS'AS", where a is such that a ■ a = 1. 
The equations of motion (I4.10p then reduce to 

C(S-BS)a + D[(S-S)l-S®S]B-ia = pa, a = aS' A S", P = 0. (6.7) 

The dot product of fl6.7P i by S = S' — iS" yields 

(S • S)(a-B-^S) - (S ■ S)(a-B-^S) = 0, (6.8) 

or, separating real and imaginary parts, 

/ (S"-S")(a-B-iS')-(S'-S")(a-B-iS") = 0, 

\ (S' • S")(a-B-iS') - (S' ■ S')(a-B-iS") = 0. ^^'^^ 

Using (S' ■ S')(S" ■ S") - (S' ■ S")^ = a^^a ■ a ^ 0, we see that this homo- 
geneous linear system of two equations admits only trivial solutions, that is 
a-B-^S' = a-B-^S" = 0. Thus a is parallel to B'^S' A B'^S", so that for 

some /3, 

a = /3(B-^S' A B-^S") = /3B(S' A S") = a'^/^Ba. (6.10) 

This last equality means that a is an eigenvector of B. Therefore, provided 
S is not isotropic, linearly-polarized inhomogeneous plane waves of com- 
plex exponential type can propagate in a deformed Mooney-Rivlin material 
only when polarized in a principal direction of the primary homogeneous 
deformation. Examples of such waves are given in section 16.41 When S (or 
equivalently, C) is isotropic, the inhomogeneous plane wave can be linearly- 
polarized in nonprincipal directions (such as in example (15.60 when a\ = 0). 
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6.3 Circularly-polarized waves, C • C 7^ 

Here, we establish the condition for the propagation of circularly-polarized 
inhomogeneous plane waves with a nonisotropic slowness bivector. From 
(14 ■2p and f l6.ip . the equations governing the propagation of such waves are 

C(C-BC) A + D{C ■ C)l-i A - D{A ■ M-^C)C = pN'^A, 

A-C = 0, A-A = 0. ^ ' 

Taking the dot product of fl6.1ip i by A yields A-B^^A = 0, or, using (13 ■4p . 

A-n^ = 0. (6.12) 

Therefore, A is orthogonal to n''^. In other words, the planes of circular 
polarization for inhomogeneous plane waves with a nonisotropic slowness 
bivector are the planes of central circular section of the B~^-ellipsoid. This 
is the same situation as for homogeneous plane waves [?]. 

From (16.111) 9 and (16.121) it follows that the projection of the ellipse of C 
onto a plane of central circular section of the B~^-ellipsoid is a circle. Note 
that there is an infinity of such C. 

We consider the case where A-n^ = and denote the corresponding 
isotropic amplitude bivector by Aq. Noting that n+ is defined by (13.31) . it 
follows that A^ is given by O §2.2] 

A^ = ae'\j A n+ + ij), a, 9 real. (6.13) 

Now, because Ag, Aq, and n"*" are linearly independent bivectors, the 
corresponding propagation bivector Cq (say) can be written as Cq = AAq + 
/xAq + i^n+. From (I6.IIP 9 we have Cq ■ Aq = /iAq ■ Aq = 0, and therefore 
/X = 0. Then from (14.31) . we have Cq.Cq = — 1 = z/^, and therefore, v is 
real and given by z/ = \/ni? — 1. Also, Cq ■ Cq = m? + 1 = 2\\a^ + i/^, and 
so |Aq;| = 1. We conclude that C can be written as 

C^ = e'\] A n+ + ij) + - ln+. (6.14) 

Using (I6.13P and (I6.14p . we compute the eigenvalue p{Nq)~'^ (say) from 
(I6.1ip i and find that it is given by 

p(iV^)-2 = A^[CA^ + D^^^F^l (6.15) 

where 

= \fV^^^^ + e^'^{\^^ - - Xf). (6.16) 
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Example 2: Circularly-polarized waves 

As an example of circularly-polarized solutions, we consider the case where 
9 = 0. In that case, we deduce from fl6.13p - fl6.16[) the following solution to 
the equations of motion f l3.8p . 

Ui = a sin0e~'^^©^ cosci;(A^Qm • x — t), 

U2 = -ae-^^Qysmu{N^m-x-t), (6.17) 
Ms = — acos(/)e~'^^o^cosa;(A^Qm ■ X — t) and p* = 0. 

Here a, u and m are arbitrary, cos(f) and sin0 are defined by (13.31) . Nq is 
given by (I6.15I) - (I6.16I) with 6 = 0, and m is the vector defined by 

m = jAn^+y/m'^ — In^ = (sin (f)+Vm'^ — 1 cos 0)1— (cos0— V^ri^ — 1 sin (f))'k. 

Note that is real, and therefore the planes of constant phase are orthog- 
onal to the planes of constant amplitude (see (14.51) ). 

The wave described by (I6.17P is circularly-polarized in the plane orthog- 
onal to the acoustic axis n+, as would be the case for a homogeneous plane 
wave. The wave is attenuated in the y-direction and travels with speed 
{uiNq)'^ in the direction of m, which is not along an acoustic axis, in con- 
trast to homogeneous waves. 

6.4 Elliptically-polarized waves, C • C 7^ 

Now we consider elliptically-polarized waves. In order to find the planes 
of polarization, we follow a procedure introduced by Boulanger and Hayes 
[7] for the propagation of finite-amplitude homogeneous plane waves in a 
deformed Mooney-Rivlin material. However, their method, which dealt with 
real vectors, is generalized here to the case of bivectors. 

We know from section 16.31 that inhomogeneous plane waves (with a non- 
isotropic S) are not circularly-polarized when the projection of the ellipse of 
C upon the plane orthogonal to is not a circle. Choosing a and b, two 
unit vectors such that (n^, a, b) form an orthogonal triad, we decompose C 
as C = (C ■ a)a + (C • b)b + (C ■ n+)n+, and see that we must have 

(C-a)2 + (C-b)V0, or C ■ C - (n+ ■ C) V 0. (6.18) 

This last inequality is equivalent to 

n±-nn± ^ 0. (6.19) 

Within this context, we seek the nonisotropic eigenbivectors A± (say) of 
the tensor n[C(C-BC)l + D{C ■ C)M-^]U such that A± ■ C = 0. We note 
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that with the Hamilton cychc decomposition (13 ■4p of the tensor B ^, we may 
write the tensor IIB"-'^!! as 

nB-^n = A^^n - ^a-^ - \-^){Un+ ® nn- + nn- ® nn+). (6.20) 

Using (16.191) . we construct the bivectors H^, given by 

H± = ^S^, H±-H± = l. (6.21) 
Vn±nn± ^ ^ 

Note that nH=^ = H=^, so that we may write flOU]) of the tensor IIB^^n as 



nB-^n = X^'^U-l{Xf-X^^)^{n+-Un+){n--Un-)[Il+(^Il- + Il-®Il+]. 

From the definition (13. 3p of and (13. 4p . we note that n"*" A n~ (parallel 
to j), n+ + n~ (parallel to i), and n+ — (parallel to k) are orthogonal 
eigenvectors of B~^. Similarly, it can be checked that the eigenbivectors of the 
symmetric operator IIB"^!! are C (parallel to H+ AH~) and the orthogonal 
bivectors A± defined by 

= ± H = -S= ± (6.22) 
Vn+nn+ vn Iln 

Now, because 11 A± = A±, we can find the amplitude bivectors A, solutions 
to (16. 2p . and orthogonal to the propagation bivector C. They are the two 
bivectors A± given by (I6.22p and their corresponding eigenvalues pN^^ are 
given by (16. 6p . or in this case, by 

pN^' = C(C-BC) + ^(C-C){(A3-2 + Ar^) 

+ (As^ - Ar^)[ ^"^"c.c'"^^ T v/K-nn+)(n-.nn-)]}. 

We note that 



p{N-^~NZ^) = -D{C ■ C){\f - Ar') V(n+-nn+)(n--nn-), (6.23) 
= -D{C-C){\f-\-,' 



fn+ ■ C) 



c c 



(n-CP 

c c 



When the bivector C and the complex quantities A^^^ are replaced by n 
(direction of propagation) and ci, C2 (speeds) for homogeneous waves, this 
equation reduces to p(c2 — c\) = D{X^^ — Aj"^) sin(/)+ sin 0~ [7], where 0^ are 
the angles between n and n''^. This in turn is reminiscent of the 'law of the 
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product of the two sines' {la loi du produit des deux sinus), for the propa- 
gation of hght through a biaxial crystal in classical linear optics, established 
empirically by Biot [20] in 1818 and theoretically by Fresnel [21] in 1821. 

From (16.231) and fl6.19p . we see that here, the eigenvalues are distinct 
and therefore neither A_|_ nor A_ is isotropic [121 §3.2]. This means that, 
as expected, the waves corresponding to these amplitudes are elliptically- 
polarized [TH] (and not circularly-polarized). In this connection, it may be 
recalled that [15] finite-a.mplitnde inhomogeneous plane waves of complex ex- 
ponential type can only be /mear/y-polarized (and not elliptically-polarized), 
when they propagate in an incompressible elastic material. Here this restric- 
tion is removed because we are dealing with sma//-amplitude inhomogeneous 
plane waves. 

Finally, using (16. ip and the fact that A± ■ A± 7^ 0, we can write pN^'^ in 
terms of C and A± as 

pN7^ = C(C-BC) + D(C ■ C) " ^ ■ (6.24) 

For homogeneous plane waves, these expressions reduce to results established 
by Boulanger and Hayes ^: replacing N^'^ with their squared wave speeds 
Ci, C2, and the bivectors C and A± by their real unit vectors n, a, and b, 
respectively, the expressions are transformed into pcf = C(n-Bn)+D(a-B^^a) 
and p4 = C(n-Bn) + D{h-M-^h). 

We now write down explicit examples of elliptically-polarized waves. 

Example 3: Elliptically-polarized waves 

In this example, we present waves propagating in a principal direction with 
attenuation in another principal direction. It is seen that for one solution, 
the wave is elliptically-polarized in the plane of the slowness bivector, while 
for the other solution, the wave is linearly-polarized in the direction normal 
to the plane of the slowness bivector. 
We take C = mi -|- ij with m > 1. 

For waves corresponding to the eigenvalue pN^'^, we have for the ampli- 
tude: A+ = a{i + imj), with eigenvalue: pN^"^ = C(m^A^ — A2) + -D(m^A2"^ — 
A^^), so that the corresponding solution is given, for real, by 

Ml = ae~'^^+^ cos ijj{mN+x — t), 

U2 = —ame~'^^+^ smu{mN+x — t), (6.25) 
Us = and p* = 0, 

where u, a and m are arbitrary (m > 1), and 

iV+ = ^p/[C{m^Xl - XI) + D{m^\-^' - \-,% (6.26) 
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For waves corresponding to the eigenvalue pNz"^, we find the amphtude 
to be: A_ = ak, with eigenvalue: pNZ^ = C{m?X\ - Xf) + D{m? - 1)A3 ^, 
and the corresponding solution is given, for real, by 

Mi = 0, M2 = 0, u-i = ae-'^^-y COS uj{mN_x-t), = 0, (6.27) 

where w, a and m are arbitrary (m > 1), and 

iV_ = ^Jp/[C{m^\l - \l) + D{m^ - (6.28) 

The waves described by (16.251) and (16.271) propagate in the x-direction 
with speed {mN±)~^, where m is prescribed and N± is given by (I6.26P and 
(I6.28p . respectively. They are attenuated in the y-direction. The wave given 
by (16.251) is elliptically-polarized in the a:?/-plane (two conjugate radii of the 
ellipse are a\ and maj); the wave given by (I6.27P is linearly-polarized in the 
z-direction. 

Example 4: Elliptically-polarized waves 

Our next example is a wave solution propagating along an acoustic axis with 
attenuation in the direction orthogonal to the acoustic axes. 

Let C = A/2n^ + y = a/2 cos 01 + + -\/2sin(^k, where 20 is the angle 
between the acoustic axes, defined by (13. 3p . After some calculation, we write 
another explicit solution to the incremental equations of motion (13. 8 p as 

= [_i ± l- 2cos2 ^ cos0e-"^±^^sina;(%/2Ar±n+ ■ x - t), 

V cos 40 

U2 = -V2[l ± ^£|^]e— cosa;(V2A±n+ ■ x - t), (6.29) 
V cos 40 

„3 = [_i ^ 1 + 2cos2 j ^.^ ^ g-^iv±y t^(v^A/-^n+ ■ x - t) and p* = 0, 

V cos 40 

where u is arbitrary, and n+ are defined by (13. 3p . and 

pNg^ = C{2X^^-Xl) + ^{Xf+Xf+{Xf-Xf)[2cos2(f)TVc^sA^]}. (6.30) 

The waves described by (16.290 propagate in the direction of the acous- 
tic axis n+ with speed {\/2N±)~^, where N± is given by (I6.30p . They are 
attenuated in the y-direction, and are elliptically-polarized. 
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Example 5: Elliptically-polarized waves 

Our last examples are inhomogeneous plane waves for which the planes of 
constant phase are not orthogonal to the planes of constant amplitude. 
Let C = mp + iq (m > 1), where 

p = cos 6^1 + sin 6'k, q = sin 6^1 - cos 6'k, < 6* < 7r/2. (6.31) 

The corresponding slowness bivector S = A^C has real and imaginary parts 
given by 

S' = mN'p - N"ci, S" = mN"p + A^'q. (6.32) 
It is found from fl6.22p that the amplitude bivectors A± are given by 

A+ = Q;(mq — zp), A_ = aj. (6.33) 

Corresponding to A_|_, we have the inhomogeneous wave, 

u = ae-'^^"-^[mqcosa;(S' ■ x - t) + psinw(S' ■ x - t)]. (6.34) 

The directions of propagation and attenuation are those of S' and S", given 
by fl6.32p . where N' and N" are the real and imaginary parts of A^, which is 
given from 06.241) by 

TTl^ — 1 

pAr-2 = C[m\p-Mp) - (q-lq)] + D —^-^[m\q-M-^ q) - (p-l"^p)] 

+ 2^m[C(p-Bq) - D— -(p-B-^q)]. 

+ 1 

This wave propagates in the direction of S' with speed |S'|^^, is attenuated 
in the direction of S", and is elliptically-polarized (two conjugate radii of the 
ellipse are mp and q). 

Corresponding to A_, we have the wave 

u = aje-"^"-^cosa;(S'-x-t). (6.35) 

Here, the respective directions of propagation and attenuation are those of 
S' and S", given by (14. 4p . where A^' and A^" are the real and imaginary parts 
of A^, which is given from (16.241) by 

= C[m^{p-Mp) - (q-Bq)] + D{m'^ - + 2im(p-Bq). (6.36) 

This wave propagates in the direction of S' with speed |S'|^^, is attenuated in 
the direction of S", and is linearly-polarized in the ^/-direction. In this connec- 
tion, it may be noted that /inite-amplitude linearly-polarized inhomogeneous 
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plane waves of complex exponential type can propagate in a homogeneously 
deformed Mooney-Rivlin material only when the planes of constant phase 
(orthogonal to S') are orthogonal with the planes of constant amplitude (or- 
thogonal to S") [n]. Here, from ([6321), we have S' ■ S" = (m^ - 1)A^'A^" ^ 0, 
and these planes are not orthogonal to each other. This is possible because 
we are dealing with sma//-amplitude inhomogeneous waves only. 

7 Summary and concluding remarks 

We have found all possible inhomogeneous small-amplitude motions of com- 
plex exponential type that can be superimposed upon the large homogeneous 
static deformation of a Mooney-Rivlin material. The results are summarized 
in Tables 1 and 2, where the solutions of exponential sinusoidal form are listed 
for inhomogeneous plane waves with an isotropic and a nonisotropic bivector 
C, respectively. In the last column of the tables, it is recalled whether or not 
each motion is possible with a finite amplitude, and the corresponding refer- 
ence is given. In the second column, the equation numbers of corresponding 
examples in this paper are given. Thus for instance, line 1.3 reads as follows: 
a small-amplitude circularly-polarized inhomogeneous plane wave of com- 
plex exponential type with isotropic propagation bivector C can propagate 
in a homogeneously deformed Mooney-Rivlin material; for this wave, the 
amplitude bivector A is isotropic and orthogonal to C; an example of such 
solution is given by equation (15. 9p : also note that A is parallel to C and that 
the complex scalar slowness is arbitrary; finally, such a finite-amplitude 
wave solution is not possible in a homogeneously deformed Mooney-Rivlin 
material (see [15J). 

Through the use of bivectors, all inhomogeneous solutions of complex 
exponential type for the problem of small deformations superposed on a large 
static triaxial strain of a Mooney-Rivlin incompressible hyperelastic material 
were obtained systematically. A great diversity and richness of solutions was 
uncovered, using the Directional Ellipse method [T6] . 

For inhomogeneous waves with an isotropic slowness bivector (^, any 
direction of polarization or plane of polarization is permitted for linear or 
elliptical polarization, respectively; also, for circular polarization the complex 
scalar slowness can be arbitrarily prescribed. 

For inhomogeneous waves with a nonisotropic slowness bivector (^, the 
secular equation was established and solved ( §6.ip . generalizing the secular 
equation for homogeneous waves [8j. Then it was seen that such linearly- 
polarized waves ( §6.2p can only propagate along one of the principal axes 
of the finite static stretch; that circularly-polarized waves ( §6.30 can only 
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Polarization Type 


Displacement 


Notes 


Finite case? 


1.1 Linear 










A-C = 0, AAA = 0. 
Example: fl5.6j) with = 


A = a2C AC, 


yes [llj 


1.2 Elliptical 


{Ae'^(^*^-''-*) + C.C.}, 








A ■ C = 0, A ■ A 7^ 0. 
Example: (15. 6p with ^ 


A = aiC + a2C A C 


no P] 


1.3 Circular 


{Ae''^(^^-"-*) + C.C.}, 

A ■ C = 0, A ■ A = 0. 
Example: (|5.9|) 


A = ttiC, 
is arbitrary. 





Table 1: Inhomogeneous plane waves, C isotropic (C ■ C = 0) 



Polarization Type 


Displacement 


Notes 


Finite case? 


2.1 Linear 




A along 


only when 




A-C = 0, AAA = 0. 
Examples: to.27p. m.m 


principal axis 
of basic strain 


S' ■ S" = 0, 
(S = NC) lllj 


2.2 Elliptical 


{Ae'^(^<^-''-*) + C.C.}, 








A-C = 0, A- A ^0. 
Examples: 

fl6.25|). fl6.29|). (16.341) 


(n± ■ C) V C ■ C 


no [15] 


2.3 Circular 


{Ae'"(^*^-"-*) + C.C.}, 

A ■ C = 0, A ■ A = 0. 
Example: (16.171) 


A ■ n± = 





Table 2: Inhomogeneous plane waves, C nonisotropic (C • C 7^ 0) 
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be polarized in one of the two planes of central circular section of the B~^- 
ellipsoid, but can propagate in directions other than the directions of the 
acoustic axes of this ellipsoid; and that all elliptically-polarized wave solutions 
( §6.41) can be obtained by generalizing a method introduced in [7] to the 
consideration of bivectors. 

Also, in contrast to the case of finite- amplitude inhomogeneous plane 
waves [m [15], it was seen (§3 and §6.41) that elliptical polarization is possi- 
ble, and that for waves linearly-polarized in a principal direction ( §6.21 and 
Example 5), the planes of constant phase need not be orthogonal to the 
planes of constant amplitude. 

Finally, it is of interest to note that the results established here can easily 
be specialized to the case of a biaxial pure homogeneous static prestrain, 
simply by taking Ai = A2 = A (say), and A3 = /i = A~^ (say). In that case, 
the tensor ni^^n may be written as [7j 

UM-^U = \-^U- {\-^ - H''^)Uk(g)Uk. (7.1) 

It follows that the computation of the amplitude bivectors and correspond- 
ing complex scalar slownesses is greatly simplified, and that the results of 
Boulanger and Hayes [Hj are directly recovered from the analysis presented 
here. 
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